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We present an exact, analytic solution of the spin dependent quantum transport problem with 
spin-orbit interaction in a one-dimensional mesoscopic ring with one input and two output leads. 
We demonstrate that for appropriate parameters spatial interference in the ring leads to a behavior 
analogous to that of the Stern-Gerlach apparatus: different spin polarizations can be achieved in 
the two output channels from an originally totally unpolarized incoming spin state. It is shown 
\Q ' that this requires an appropriate interference of states that carry oppositely directed currents. We 

find that spin polarization is possible for several geometries, including the case when the device 
is not symmetric with respect to the incoming lead. A clear connection is established between 
the Stern-Gerlach like property of the device and the relevant Aharonov-Casher phases in the loop 
q ' geometry. 
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I. INTRODUCTION 



Electronics based on the spin of the electron represents a new direction of development (spintronics) [1]. In order 
to utilize this additional degree of freedom as an (either classical or possibly quantum) computational resource it is 
required to develop a controllable way of manipulating spins. One of the most promising mechanisms that can be used 
C*~) ■ for this purpose is the spin-orbit interaction [2J in semiconductor materials, which can be controlled by an external 
electric field For semiconductor nanostructures, where the mean free path of the electron can be much larger 
than the size of the device, quantum mechanical interference can lead to a new class of spin-sensitive devices, such 
as quantum gates [ij . On the other hand, even if a whole spintronic apparatus does not use the quantum mechanical 
\q , nature of the electron for information processing, some parts of it still may rely on interference phenomena, similarly 

1 to the polarizing device discussed in this paper. 
""^ | Quantum rings @, 0, HJ or loops 0] (i.e., ring shaped objects where quantum mechanical interference plays 
an important role) made of a semiconductor material, have been shown to have remarkable spin transformation 
I ; properties @, fl i |, M HI HI El M El ■ 

This is partially due to the geometry of these devices, as the incoming 
electrons are forced to split into two different spatial parts that interfere at the output, while the spin-sensitive 
interaction introduces an additional effect to be taken into account. Consequences of the interplay between spatial 
interference and the so-called Rashba-type [2| spin-orbit interaction was recently observed in HgTe nanorings, where 
an external magnetic field was also present |l6l ] . More generally, the spin degree of freedom in quantum interference 
[TtI ITU ITjjl I20I can play a prominent role in the development of a spintronic network based on various spin-sensitive 
devices [H M SI ■ 

Let us recall that already Bohr and Mott pointed out [26|, that in contrast to atoms, one can not spin-polarize 
electrons in an inhomogeneous magnetic field. We consider here a three-terminal quantum ring, where electrons 
entering in a totally unpolarized spin state become polarized at the outputs with different spin directions. This device 
can be deemed in a certain sense a spintronic analogue of the Stern-Gerlach apparatus 25]. A related polarizing effect 
was recently predicted in a Y-shaped conductor which was a consequence of scattering on impurities [27| . Note that 
this is a very different physical mechanism from the coherent spin transfer to be discussed here. Our model is based 
on an exact, analytic solution of the spin dependent transport problem and thus provides a clear physical picture 
of a process where fundamental polarization effects as well as nontrivial spatial-spin correlations, entanglement or 
intertwining (28j appear (29j . In addition, our treatment allows us to determine the parameter values for which the 
device is rcflcctionless, i.e. perfect polarization at the outputs takes place without losses. 

In the present paper we demonstrate that the physical origin of the polarizing effect we found earlier in Ref. |25| is 
essentially spatial interference: At a certain output junction the spatial parts of one of the eigenspinors representing 
clockwise and anticlockwise directed currents interfere destructively, leading to the transmission of the other orthogonal 
eigenspinor in the output lead. This effect is visualized by plotting the spatial dependence of the spin direction along 
the ring. We show that there are several, not necessarily symmetric positions where such destructive interference 
takes place. All of these configurations possess the polarizing_property. We also determine the conditions for the spin 
polarization effect in terms of the Aharonov-Casher phases [30] gained by the two orthogonal eigenspinors. 
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II. THE MODEL OF SPIN DEPENDENT SCATTERING IN A RING 



The one-dimensional Hamiltonian of an electron moving on a ring situated in the x 
Rashba spin-orbit interaction is given by[ll], [3l| 



y plane in the presence of 



h = nn 



^ + 2T2 ( ^ C ° S ^ + ^ Sin ^) ~h 



(1) 



where tp is the azimuthal angle of a point on the ring, h£l — ft 2 /2m*a 2 is the dimensionless kinetic energy of the 

electron, a is the radius of the ring, m* denotes the effective mass of the electron, and u> = ehE z j {\p2m*ac} is the 
frequency associated with the spin-orbit interaction, with E z being a static electric field perpendicular to the surface 
of the ring. The parameter uj can be tuned with an external gate voltage 0. Apart fr om constants, this Hamiltonian 
([T]) is the square of the sum of the z component of the orbital angular momentum operator L z = — and of 
fjSV, where S r = o>/2 is the radial component of the spin (both measured in units of Ti). H commutes with the z 
component of the total angular momentum K = L z + S z , as well as with Sg v = S x sin 9 cos if + S y sin 9 sin if + S z cos 9, 
the spin component in the direction determined by the angles 9, and <p, where 9 is given by tan# = — ui/fl. Therefore 
a basis can be constructed which are simultaneous eigenfunctions of the operators H, K, and Se v - In the |t), ||) 
eigenbasis of S z we can find these in the form Q: 



(2) 



obeying 



The energy eigenvalues are 



Kip (k, ip) = kiJj («, tp) , 
Se 9 ip(K,(p) = ±-i/)(k,<p). 



= HQ, 



1 $(m) 



K ~2~ 



2n 



= mi [k 2 + (-iy 



U3 



1/4] 



(3) 



where fj, = 1, 2, w = \/l + u> 2 /Q 2 and $^ = — ir [1 + (— 1) M w] is the Aharonov-Casher phase [30]. In a ring connected 
to leads the energy is a continuous variable - since (ft ± 1/2) is no longer an integer as it is in the case of a closed 
ring - and the possible values of k are the solutions of equation (J3]) for E^ = E: 



K J = ~ 



Z Z7T 

= (-ir +1 (w/2 + (-iy q ), 



(4) 



where fi,j — 1,2 and q — y (w/217) 2 + E /Ml. The energy eigenvalues of the Hamiltonian are fourfold degenerate, 
thus the state of the electron in the different segments of the ring (see figure [1]) for a given E is a linear combination 
of the corresponding four "4>{k^, (p) eigenstates 



%(cp)= i = i,n,in, 

where the ratio of the components of the spinors in ([2]) is given by 



(5) 



tan6» ( ^/2- - (1 + {-if w) . 



Since tan0^ 1 )/2 = — cot#( 2 )/2, we can express the two eigenspinors with = 9: 



(6) 



,i,,J lK >fe-^/ 2 cosf 
e l(fi / 2 sin f 



(7) 
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FIG. 1: The geometry of the device and the relevant wave functions in the different domains. The parameter <p is measured 
from junction 3 in counterclockwise direction. 



The stationary states of the complete problem (ring and leads) , can be determined by imposing continuity of the 
wave functions at the boundary of the different domains. Using local coordinates as shown in figure [TJ the incoming 
wave, ^(a^), and the outgoing waves ^(xi), ^'2(^2) are built up of linear combinations of spinors with spatial 
dependence e ikx etc. corresponding to E = h 2 k 2 /2m*: 

= ($)e <te », (9) 

where n = 1,2. Our aim is to determine the transmission properties of the ring, i.e. to find the elements of the 
transmission matrices, which arc defined in the following way 

r< "> $)-(?)■ ,10) 

where n = 1, 2 label the two outgoing leads. 

In order to obtain the transmission matrices, the 12 coefficients have to be determined. This can be done via 
requiring continuity of the wave functions, and vanishing net spin current densities (Griffith conditions) pi. Ull [32l [33j 
at the three junctions. 
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According to the detailed calculations presented in the Appendix, the elements of the transmission matrices are: 



T, 



(n) 
11 



-e • 2 



» 



» 



T 



(n) _ 

n ~ 

(n) _ 
IT 



' 2 sin - cos - 
V 2 2 



(fci n) +/4 n) )-(fti n) -^ n) )' 



e J n , 



T 



(n) 



u 



a ,22, 
-e 2 



sin 



+ /4 n) 



cos 



(™) l(") 



where 



h{' = -kae 1 ^ 11 sin(g(27r - 72)) sin (q (72 - 71)) , 
h{ 2) = fcae l f^e- 4t, " r sin(( 7 7i)sin(q(72 -71)), 
4 n) = iqe^ ( e -" U7r sin(g7„)-sin(g(2 7 r- 7 „))) , 

and y is given by (|A53|) . 

The reflection matrix R is found to be diagonal in the {||), ||)} basis of S z : 



R 



11 



R 



11 



R\\ = - — {— iq 2 sin (2qTr) + ik 2 a 2 sin (q (72 — 71)) sin (q (2tt — 72)) sin ((771) 
.V 

-gfca [sin (q (2ir - 71)) sin (qji) + sin (g (27r - 72)) sin (qj 2 )}} - 1, 
i?XT = 0. 



(11) 



This matrix describes the losses in the efficiency of the spin transformation, as the sum of the norms of the outgoing 
and reflected waves should be equal to the norm of the input. 



III. ANALYSIS AND VISUALIZATION 



When the incoming electron is not perfectly spin-polarized, i.e. its spin state is a mixture, which - instead of a two 
component spinor - should be described by a 2 x 2 density matrix g ini then we can readily generalize equation (|10[) 
to obtain 

where g^- 1 ' and g^ are the output density matrices in the respective leads. 

Considering a completely unpolarized input, i.e. gi n being proportional to the 2x2 identity matrix, in order to 
get polarized outputs, the relevant output density operators should be projectors (apart from the possible reflective 
losses): 

1 T («) (t^ =r) n \<l> n )(<t> n \- (12) 

The non-negative numbers r\\ and 772 measure the efficiency of the polarizing device, i.e. 771 + 772 = 1 means a 

reflectionless process. Equation (p~2|) is equivalent to require that the determinants of (T^ n A vanish. According 
to equations pip there are two different conditions for each transmission matrix to satisfy this requirement: 

h[ n) ± /4" } = 0, (13) 

where n = 1, 2 indicates the two output juctions. It can be shown that only the following two cases lead to nonzero 
transmission at both outputs: 

+ h? = 0, 

h (2) ~h {2) = 0, (14a) 
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tf> - = 0, 

h^+h^ = 0. (14b) 
After substitution we obtain two equations for cos (wit) and sin (ww) in both cases: 

cos (w) = Sin(g(2 , ?r 7 l)) = . ^ (15a) 

sm(Q7i) sm(g(27r-7 2 )) 

. fca sin (g (2?r - 72)) sin (g ( 72 -71)) fca sin (571) sin (g (72 - 71)) 

sm(w7r) = ± — = ± . -, , n rr , (15b) 

q sin (971) q sm {q (2-7T - 72)) 

or with the Aharonov-Casher phase 

cos$ (.) = _sin(g(27r~7i)) = sm( TO ) 

sm(g7i) sin (9 [2tt - 72)) 

sin $M _ ± , r y ka sin (g ( 2?r ~ 72)) sin (g ( 72 - 7l )) 
9 sin (571) 

= ± (-1)" ^^( g 7i)sin( g ( 7 2-7i)) ; (16b) 
g sm (g (27r - 72)) 

where the plus sign in (|15bj) corresponds to equation ()14al) . while the minus sign is for the case given by equation 
(|14b|) . From (|15ap and (] 1 5b[) for both signs we find 

sin(g7i) = ± sin (g (2ir - 72)) , 
sin(g7 2 ) = ± sin (g (2ir - 71)) . 

The solutions of these equations for a fixed 72 are 

71 = 27r - 72 ± mn/q, (17) 

or for a fixed 71 are 

72 = 2tt - 71 ± frr/g, (18) 

where m, I are nonnegative integers which ensure 0< 71,72 <2tt and 72>7i- The m,l = cases correspond to a ring 
the outgoing leads of which are symmetric with respect to the incoming lead. It was demonstrated that for such a 
ring one can find lines in the {72, w/f2, ka} space along which the conditions (|15ap and (|15b[) can be satisfied [25| . i.e. 
the ring polarizes a completely unpolarized input. Polarization occurs with equal 771 = r/ 2 = rj/2 transmission in both 
outputs. Parameter combinations, for which the transmission probability r\ is unity can also be found. 

From equations (|17[) and (|18[) we see that we can extend the polarizing property to asymmetric geometries. The 
asymmetric positions of the two output leads for which the condition for complete spin polarization is satisfied 
are ±lir/q and ±mir/q angles away from the symmetric ones. For proper combinations of the uj/VL and ka pa- 
rameters the asymmetric ring also produces polarized outputs with equal transmission probabilities 771 = 772 = 



128g 2 fc 2 a 2 



/ |y| 2 . We note that this is an important generalization of the results of Ref. [25| . There are several 

appropriate positions for the output leads, the symmetric case is just one of them. 

The output spinors are the eigenstates \<j) n ) of the transmitted density matrices iT' n ' (T^) , which correspond to 
the nonzero eigenvalues given by T) n . Focusing on the case of equations (|14a|) . these eigenstates read 



ii\ / e " 2 sm S \ 1 ,2\ / e 1 2 cos % \ , 1 _ N 



cos 2/ \ e 2 sm 2 

These results describe the connection between the strength of the spin-orbit coupling (encoded in 9), the geometry of 
the device and its polarizing directions. Note that these spinors are in general non-orthogonal, their overlap is given 
by (0 2 | 1 ) = i sin 9 sin (72 — 7i)/2. For the other case given by equations (|14b[) we have: 

m = ( e lS co i), \4)=( e ~Z s H)- (20) 

1 ' \ e'i sin? i 1 ' V— e l 2 cos S / 
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FIG. 2: Transmission probability r\ = 2r\\ — 27/2 at the two outputs of an asymmetric ring as a function of ka for w/fi = 2.27. 
This figure corresponds to 72 = 47r/3 and 71 = 2-7T — 72 — 6n/q. The dots mark the points where perfect polarization occurs. 

We can see from (fl9|) and (f20|) that the output spin states in both cases are the two eigenspinors of the Hamiltonian 
at the positions of the output junctions. For a given output lead, the two eigenspinors are interchanged in the two 
cases. 

Figure [2] shows this transmission probability as a function of ka for u/Q = 2.27. The dots on the curve mark 
the points where perfect polarization occurs. Since the angle 71 is a function of ka, the dots correspond to different 
asymmetric configurations. It can be seen that even for an asymmetric ring, with appropriate parameter values, 
complete output spin polarization can be achieved with practically zero reflective loss. 

Now we investigate the physical origin of this polarizing effect. To this end we consider a completely unpolarized 
input taken as the following equal weight sum of pure state projectors 

ftn = \ (M„) (i>in\ + K 2 „) <V>in|) ■ 

Here ipi n = V'i'(O) = 1P2 (0) (l 1 — 1)2), given by equations jJJ and ([8]), are the eigenspinors of the Hamiltonian at the 
position of the incoming lead (3). The density operator in the different sections of the ring is then 

ft(^ = i(|*-(^)<^M| + |^M)<* 4 2 (^)|) i = 1,11,111, (21) 

where 

J'=l,2 

3=1,2 

are the spinor valued wave functions of the electron in the different domains of the ring for the pure inputs ip} n an d 
ipf n respectively, with 

N ? O) = < e "^ M=l,2. 

J'=l,2 
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We can see that for these inputs the wave functions in the ring contain only two of the four eigenstates of the 
Hamiltonian, those which have the same spinor part. 

By calculating the spin current densities corresponding to the "J^ ((f) states appearing in (f22|) we obtain 



2^ 



-1 



r (c 



)]=(-!) 



(23) 



By examining (|23p we find that (<p) and ^ 2 (<£>) represent oppositely directed (clockwise and anticlockwise) spin 
currents in each section (identified by the index i) of the ring, since Jf^ and J^ 2 have opposite signs. The overall spin 
current densities - containing both clockwise and anticlockwise directed currents - which correspond to the input ip^ n 



J? 



2g(-ir(|<f 

2 q (-ir(\< 1 \ 



a: 



s2 



-2Re 



COS! 



(2 



^12 ■ 



(24) 



We note that the disappearance of the cross terms in ([24]) is due to the fact that tan (9 = —ui/Cl. 

The output spinors given by (Til?]) and (f2T)]) suggest, that in order to obtain a polarized (pure) state at a given 
output, we need one of the one-dimensional projectors of (|21[) to vanish, and the other one to remain nonzero at that 
point of the ring. In order to have different polarized spin states in both outputs, the two projectors need to vanish 
at the different output junctions. One of the possible ways to achieve this is to have \P} (71) and (72) being zero, 
which happens if the spatial parts of these wave functions are zero (see equations (1221) ): 

Nj( 7l ) = 0, Nfj( l2 )=0, (25) 

indicating destructive interference at the given output. By exchanging /i = 1 and /i = 2, we can describe the other case 
of polarization. Condition (|25[) can be satisfied if 

Kl| = K,2|> | a //,l| = |«?Z,2 1 . ( 26 ) 

which, by using equation (|A52[) . can be shown to be equivalent to equations (|14a[) . (Exchanging /i = l and /i = 2 in 
(f2"5| leads to (|14b[l ). Equation (j2"6")l implies that the spin currents J} and Jfj given by (fM]) vanish as a consequence 
of the interference of oppositely directed currents corresponding to states of the same spinor parts. We can see that 
the requirement for spin-polarization given by (|12|) . has a very clear physical interpretation in terms of destructive 
interference and vanishing spin currents. 

Figure [3] and the corresponding animation show the stationary spin directions of the electron along the ring for a 
completely unpolarized input, for parameter values which ensure perfect polarization at the outputs (in the case given 
by (|14a[) for a symmetric ring). Light and dark (blue online) arrows correspond to ipj n and ipf n , respectively. The 
length of the arrows as well as the curves with the corresponding colour on the upper graph show the probabilities 
of finding the electron at the given point on the ring. The dashed lines in the upper graph mark the outgoing leads, 
where one of the two probabilities becomes zero, leading to the output of the other spinor, given by equation (|19[) . 
We note that spin transformation in this case is a rotation around the z-axis by an angle pertaining to the given point 
on the ring. 

It is also interesting to see how polarization is produced if we decompose the incoming perfect mixture as an equal 
weight sum of the eigenstates of S z 

ftn = ^(IT) (Tl + U) UD- 

Figure f?] and the corresponding animation show the stationary spin directions along the ring for such an input and 
for the same parameter values as in figure [3l Light and dark (blue online) arrows correspond to inputs ||) and ||), 
respectively. The outgoing arrows (green online) represent the output spinors given by equations (|19p . The spin 
direction of the electron in the two branches of the ring is illustrated on the two Bloch spheres [35| above the ring, 
where the length of the black arrow represents the purity of the given state. When the arrow reaches the surface, 
the spin state is pure, otherwise it is mixed, zero length meaning a perfect mixture. The dots on the Bloch spheres 
(green online) indicate that at the positions of the output junctions the states are pure. The animation shows that 
the |t) and ||) inputs are rotated into the same direction at the outputs of the ring, resulting in the same pure states 
as those seen in figure [3] 

The role of spin-orbit interaction in the polarizing process can also be seen in figures [3] and [4] The spin-sensitivity 
of the problem leads to a symmetry 36] of the stationary solution that is more complex than the pure geometrical 
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FIG. 3: The stationary spin directions of the electron along the (symmetric) ring for a completely unpolarized input given 
by Qin — | (\4>ln) (4>ln \ + \4>in) (4>in\) with 72 = 47r/3, id /Q = 3.05, fca = 1.38 [34[ , which ensure perfect polarization in the 
case given by equation (|14a[) . Light and dark (blue online) arrows correspond to ip] n and tpf n respectively. The length of the 
arrows as well as the curves with the corresponding colour on the upper graph show the probabilities of finding the electron at 
the given point on the ring. The dashed lines in the upper graph mark the outgoing leads, where one of the two probabilities 
becomes zero, resulting in the output of the other spinor as a pure state. The two outputs in this case are given by equations 
©. 

mirror transformation. Note that without spin-orbit interaction destructive interference for a given spin direction 
would imply that the orthogonal spin component of the wave function is also zero at that point. Consequently placing 
the output junctions in such positions would mean zero transmission probability. This can also be seen by considering 
that the rotation of the |f) and |J.) spinors shown in figure [4] is due to spin-orbit interaction, for ui — (or <f> = 0) 
their direction in the ring would not change, i.e., they could never precess into the same direction. 

Figure and 2] show that besides the actual output junctions, which are situated symmetric with respect to the 
incoming lead, there are additional points on both branches of the ring where the state of the electron is pure. These 
points are situated asymmetric with respect to the incoming lead, and can be determined from equations (|17[) and 
(jTSJ . If the output leads are put into those positions, the outcoming spins states are the ones given by (fTT)l) . 
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FIG. 4: The stationary spin directions along the ring for a completely unpolarized input given by Qi n = | (| j) (|| + \[) (||), for 
the same parameter values as for figure [3] Light and dark (blue online) arrows correspond to inputs ||) and ||) respectively. 
The outgoing arrows (green online) represent the output spinors given by equations ([19} ■ The spin state of the electron in left 
(right) branch of the ring is illustrated on the Bloch spheres [35| above the respective part of the ring. At the positions of 
the output junctions the black arrows reach the surface of the spheres - denoted by the dots (green online) - indicating that 
the state at those points is pure. The animation shows that the ||) and |J.) inputs are rotated into the same direction at the 
outputs of the ring. 



IV. CONCLUSIONS 

We considered a three-terminal quantum ring with one input and two output leads, which for appropriate parameter 
values acts as a spin polarizer, similarly to the Stern-Gerlach apparatus. We presented a detailed analytic solution of 
the spin-dependent transport problem and provided the physical interpretation of the process: For both symmetric 
and non-symmetric geometries, polarization is due to spatial interference. At a given junction this interference is 
destructive for a certain spin direction, while constructive for its orthogonal counterpart, which, consequently is 
transmitted into the output lead. 
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APPENDIX 



Considering the input junction, the continuity and Griffith conditions read: 

* 3 (0) = */ (0) = 9 m (2tt) , 

J 3 (0) - Ji (0) + J ni (2tt) = 0, 



(A27) 
(A28) 



respectively. Analogous equations can be written for the other two junctions. The appropriately normalized spin 
current densities in the leads are given by 



and in the ring by [111 ] 



Ji O) = 2Re 



(A29) 



(A30) 



where er r (tp) = cospa x + smipa y , i — I, II, III and I — 1,2,3. For the sake of definiteness, we present the details for 
the incoming junction (3). The results for the other junctions can be obtained in a similar manner. 
First, we simplify equation (|A28jl by using (|A29p and (|A30p 



0*3 



0X3 



_ 0*£ 

. 3 =o 0^ 



+ 



9* 



in 



ip=0 



dip 



0. 



Substituting the wave functions (JSJ) and @ into (|A27|) we get 



The same substitution in (|A31|) yields 



E-;t<, "1"; 



(A31) 

(A32) 
(A33) 

(A34) 
(A35) 



where f3j = (k^ 1 — |) /fca and 7^ = (k^ 1 + |) /fca. Since our aim is to determine the afj coefficients, we eliminate rj 
and rj_ 



m, 3 ) 3 u ' 



EK+« 

E [# («? 

V,3 

E[^(<, + ^^-) + < 



2/ T , 



< = 2/ X . 



/' ,3 



After substituting the m^ 1 and v% spinor components given by (|7|) and © into (|A36|) - (|A39|) we find 



J77j e^ 2,r ] cos 



E ( a L + 4 

E/ 1 I 1 IK 1 



) c °4 + ( 



Qjj + "-,„.,(" 2z ] sin ■ 



( a L" + a 



L, 3 -e^ 27r lcos. 



0, 
0, 



(A36) 
(A37) 

(A38) 
(A39) 

(A40) 
(A41) 
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COS ■ 



2 iKj27T 

a/77 j e 3 



l I,3 



sin ^ ^ = 2/t, 



"3 s g f* - 2 h- 



Notice that certain terms can be cancelled out by using simple trigonometric identities, giving: 

3 

E{ a h+ a ni^>) = 0, 



(A42) 
(A43) 

(A44) 
(A45) 



,2 „ik;27t 



[(2k] - cos 9) (a) 3 + a^j.e'^j + 2kaa 1 I>j - sm9 (aj 

3 

[( 2k j + cos0 ) (o?j + a?77je iKj227r ) + Zkaaj j - sin (9 (a^ + "]„.,' "'^ 



= 2d 1 , 
= 2d 2 , 



where 



ri 1 = 2fca ( cos -it +sin-/| 



d = 2ka ( sin -/ T - cos 



The sums in equations (|A46[) and (j A47|) can be simplified using (jA44j) and (|A45[) 

[(2«) + 2ka) a}j + 2 K y mj e iK > 

3 

[(2/t 2 + 2ka) aj tJ + 2 K 2 J aj II ^e lK > 



= 2d 1 , 
= 2d 2 . 



(A46) 
(A47) 

(A48) 
(A49) 

(A50) 
(A51) 



Thus, the equations originating from the continuity requirements at the incoming junction (3) (see figure Q} split into 
two separate systems for fi = 1,2. 

The other two junctions lead to four additional equations and consequently we have to solve six equations for six 
unknowns for each [i. We can start by expressing a\ u ■ by a\ ■ from equations (|A44[) and (|A46[) . and also by a}j ■ 
from the equations for junction 2. Using these two expressions for a] n ^ we obtain a relation between a\j j and a] ,-. 
From the equations of the first outgoing junction we can also express a}j ■ in terms of a} j. Finally, we can use the 
two different expressions for a\j j to calculate a] j . This automatically determines all the other coefficients (since we 
have these expressed by a] j)- The [i — 2 case can be solved analogously. The twelve coefficients read 



% I,3 = 



2d' 1 
V 



{kW^ 1 ^™ sin (q (<y 2 - 7l )) sin (q (2tt - 72 )) 



-iqka 



i(-i)" +3 97 1 sin ( g ( 27r _ 7l )) + e »(-i)" +J 972 gin (g ( 2 7T - 72 )) 



l U,3 



_ q 2 (^{-iy w . + e i { -ir+i q 2^y 



+ifca 

1 r 
2~q 



i(-i)"+^7 le i(-i)"- sin (g7l) + e i(-ir«<m sin (5 ( 2?r _ 72 ))j | , 
(2, + (-1)^ fca) oj&j + (-1)^' fcae 2 ^ 1 ^^ , 



(A52) 
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where 

y = ik 3 a 3 [sin (2q (tt - 72 + 71)) + sin (2q (tt - 71)) - sin (2q (tt - 72 )) - sin (2qn)] (A53) 
—2qk 2 a 2 [cos (2g (7r — 72 + 71)) + cos (2q (tt — 7l )) + cos (2g (it — 72)) — cos (2qw)] 
+4gfc 2 a 2 cos (2<j"7r) — \2iq 2 ka sin (2g7r) + 8g 3 [cos (wit) + cos (2g7r)] . 

(n) 

The elements of the transmission matrices can be determined from the continuity equations t\ — 
Sa^j.e'W"^) 7 "^ and t 2 n) = E 4/^' + ^ 7 ""i at the outgoing junctions (1 and 2) yielding (11)1 . 
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